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The lattice spacing comparator estab- 
MshcJ HC the National Institute of Stan- 
dards and Technology to measure the 
lattice spacing differences between 
nearly perfect crystals j^ described in 
detail, L,atticc pacing differences are 
inferred from the measured differences 
in Bragg angles for different crystals. 
The comparator is a two crystal spec- 
trometer used in the nearly nondisper- 
sive geometry. It has two s-ray sources, 
two detectors, and a tluvicc which per- 
mit.'; rcmtite interchange of the second 
cryilal sample. A sensitive heterodyne 
interferometer which is calibrated with 
an optical pulygon is used to measure 
the Bragg angles. The ciysiaLs are man- 
ufactured with nearly equal ihicknesses 
so that the recorded profiles cTihibit 
pcndelliisung oscillations which permit 



more precise division of the X-ray 
profiles. The difference in lattice spac- 
ing between silicon samples used at 
PbjTiikalisch-Tcchnlsche BundesarLsialt 
(PTB) and NIST has been measured 
with a relative uncertainty of t x 10'*, 
This measurement Ls consistent with ab- 
solute lattice spacing measurements 
made a1 PTB and NIST, Components 
of uncertainty associated with system- 
atic effects due to misallgr^ments arc 
derived and estimated. 
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1. Introduction 

Absolute measurements of the crystal lattice 
spacing of silicon have been reported by two labo- 
ratories [1,2,3,4] and are in progress in at least two 
other laboratories [5,6], These experiments present 
a very difficult metrological challenge and require 
extreme care in order to achieve a relative uncer- 
tainty of several parts in 10^. The output of these 
experiments is the lattice spacing of a particular 
siliom specimen in terms of an optical wavelength 
standard. 
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Crystals having well defined lattice spacings and 
whose lattice spacings are known with an uncer- 
tainty in the 0.1 to O.Ol ppm range (lppm = 10"*) 
are needed in the field of metrology (i.e., x-ray and 
gamma-ray wavelength measurements). However, 
the effort required to directly compare a large fam- 
ily of crystal samples to an optical wavelength stan- 
dard is prohibitive. Thus, several techniques and 
devices to more rapidly measure lattice spacing dif- 
ferences with a relative uncertainty below 0.1 ppm 
have been devised. This paper describes in detail 
the precision lattice spacing comparator (Ad com- 
parator) established at the National Institute of 
Standards and Technology and reports the results 
of a careful comparison of Si samples which were 
measured absolutely by the x-ray interferometer 



technique at the Physikalisch-Technische Bunde- 
sanstalt (PTB) and NIST. The measured lattice 
spacing difference has a relative uncertainty of 
1.4 X 10'* and is, at present, consistent with the ab- 
solute measurements reported by the two laborato- 
ries, 

2. Principle of the Measurement 

The principle of the NIST lattice-comparison 
measurement is illustrated by the geometrical lay- 
out of the spectrometer (Fig. 1). This geometry, 
which is similar to one suggested by Hart [7], is a 
conventional double-crystal Laue-case spectrome- 
ter. Ando et al. [8], Becker et al. [9], and Hauser- 
mann et al, [10] have also used similar geometries 
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Fig. 1. Schematic diagram of the NIST lattice spacing comparator. Si, S»— x-ray sources; Lt, Li— shutlers; Aj, A;— apertures; Ci, 
Q— crystals; Dj, Dj— deteaors; /, / ' /'— x-axis projeelions of the i-ray paths from Ai to d. Q to Ci, Q to Ai, respectively, a) Top \iew 
{x, y ptanc), b) side view (*, z plane) 
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to precisely compare lattice spacings. As will 
emerge in the discussion which follows, the NIST 
instrument has certain new features which enhance 
its generality and sensitivity. The spectrometer has 
two sources, S[ and Sj, two crystals, Ci and C? and 
two detectors, Di and D;, and is used to record 
only nearly nondispersive profiles. The x-ray beams 
are defined by precision source, Ai, and detector, 
Ai, apertures and there are computer controlled 
shutters, Li and L^, near the sources and the detec- 
tors, respectively. 

By using the spectrometer in the nearly nondis- 
persive mode, the large intrinsic line widths associ- 
ated with characteristic x-ray emission do not 
contribute appreciably to the width of the recorded 
profiles. The source shutters permit isolation of the 
two x-ray paths so that doubly diffracted profiles 
from one source can be recorded without interfer- 
ence from singly diffracted profiles from the other 
source. There are three apertures in front of each 
detector (Fig. lb), one above, one on, and one be- 
low ihe plane of dispersion. An indexed shutter se- 
lects one of these apertures depending on whether 
lattice spacing difference scans or crystal alignment 
scans are being recorded. For lattice spacing differ- 
ence scans, apertures on the plane of dispersion 
are used for both sources and the detectors. As 
crystal 1 is rotated through the reflection, the 
source shutters alternately pass and block the two 
x-ray paths, path (-••) and path ( — — ), so that the 
two x-ray profiles are essentially recorded simulta- 
neously. Thus errors due to drifts of the crj'stals or 
the angle interferometer are eliminated. For crystal 
alignment scans, only one source is used and the 
detector shutter alternately permits radiation to 
pass through the apertures above and below the 
plane of dispersion as crystal I is scanned. The an- 
gular shift of these two out-of-plane profiles is re- 
lated to the tilt misalignment of the planes of the 
first crystal with respect to the planes of the second 
crystal. By careful alignment, the contribution to 
the uncertainty due to crystal tilt misalignment can 
be made small (see Appendix A). 

General features of the comparator operation 
may be understood as follows; If the crystal planes 
near the ends of the first crystal are parallel and 
the lattice spacings of the first {di) and the second 
(d2) cr>'Sials are equal (rfi = rfi), then the profiles in 

the path (- • ■) and the path ( ) will peak at the 

same angular setting of the first crystal. However, if 
di^di, then there will be an angular offset, ^i be- 
tween the two profiles which is a measure of the 
lattice spacing difference, i^d =£f: — rfi, between the 
two crystals. If Bi and Bi are the Bragg angles of 



crystals 1 and 2, respectively, then 0i.) = 2(ft — 0i) = 
2Afi. If irf/rf is small (<iO"') and the required 
uncertainty in Ad/rf is ~ 10"^ then it is appropriate 
to use the differential form of the Bragg equation 
to express Arf/d. 



Adid = -(A0) coiB= -(^ii/2) cote 



(1) 



In our comparison scheme, the first long crystal 
serves as a temporary (a few hours) reference. Pro- 
files are recorded using one of the small second 
crystals. Then another small second crystal is 
brought into the position where the x-ray beams 
diffracted by the first crystal intersect. The first 
crystal is again scanned and profiles are recorded. 

By subtracting the angular offsets, j8, measured 
with the different second crystal samples, the prop- 
erties of the first crystal such as the absolute lattice 
spacing, the effects of variation in lattice spacing, 
and parallelism of crystal planes at the two ends of 
the crystal are eliminated. One needs only assume 
that the properties of the first crystal are constant 
during the time of measurement of the two differ- 
ent second crystal samples. The equation for the 
difference in crystal lattice spacings for the pair of 
small crystals becomes 



Ldld,= -(Ae) cotft= -0t2) cotfl,. 



(2) 



where ft = the Bragg angle for the standard sample, 
9^ = lhe, Bragg angle for the unknown sample, d^^ 
lattice spacing of standard sample, da = lattice spac- 
ing of unknown sample, Ad~dt-da, ^ = 2(ft-flt) 
= angular offset for the standard sample, 
J3u = 2(0u - flj) = angular offset for the unknown 
sample, and ^ = ft - ^^ = 2(0^ - 0^} = 2A 8. 

In Appendix A the equations for calculating Ad 
in ideal and misaligned geometries are derived. For 
small lattice parameter differences and well 
aligned x-ray beams and crj'Stals, the more involved 
equations derived in the appendix reduce to Eq. 
(2). The forms of the corrections associated with 
misalignments are explicitly given and their magni- 
tudes are estimated for our particular spectrome- 
ter. This appendix is hkely to be most interesting to 
the precision x-ray specialist. Throughout this pa- 
per and particularly in the appendix the path (- ■ ■) 

is referred to as the path (-) and the path ( ) 

is referred to as the path ( + ). 

Since X rays from the two paths interrogate the 
same area of the second crystal, small variations in 
lattice spacing along a crystal can be measured with 
this scheme. The second crystal samples reside on 
a precision slide which allows easy interchange of 
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crystal samples. A typical data sequence involves 
many changes of crystal samples and permits com- 
pensation for drifts of the crystals and the angle 
measuring spectrometer. 



3. Experimental Apparatus 

The lattice comparison spectrometer rests on a 
1 m X 2 m cast iron plate which is isolated from 
building vibrations by three passive air bags. The 
x-ray tubes, detectors, and shutters are rigidly at- 
tached to this plate. A smaller plate (61 cm x 122 
cm) on which the laser, receiver, and precision 
spectrometer are mounted rests on vibration damp- 
ing rubber feet on top of the large cast iron plate. 
The precision spectrometer is constructed on a 35 
cm X 74 cm X 2S cm thick cast iron U channel. 
The angle interferometer, rotary table, and transla- 
tion table reside on top of the inverted U channel 
and the drive arm coupled to the rotary table is 
conveniently located underneath the channel. The 
important dimensions which define the scale of the 
spectrometer are: 1) x-ray focal spot to source slit 
= 6 cm; 2) source slit to x-ray tube pivot = 25.4 cm; 
3) x-ray tube pivot to first crystal = 8.25 cm; 4) first 
crystal to second cr\-stal = 8.25 cm; and 5) second 
crystal to detector slit = 28^ cm. 

The radiation sources are identical silver or 
molybdenum x-ray tubes operated at K=40 keV 
and / = 15 mA and having a I mm x 1 mm focal 
spot (in projection). The tubes are cooled by tem- 
perature regulated water so that the spectrometer 
temperature remains constant and uniform. The x- 
ray sources pivot about the point of intersection of 
the x-ray beams which lies between the sources and 
the first crystal. Bragg angles from 10° lo 22° are 
accessible. The vertically defining slits in front of 
the sources have a 1 mm opening. The center of 
these slits and the middle slits near the detectors 
are carefully placed on the plane of dispersion 
(plane perpendicular to the axis of rotation) to 
within 0.1 mm. The x-ray tubes are positioned so 
that the focal spots fill the source and detector slits 
uniformly. The lead shutters near the source slits 
are computer controlled to block or pass the x-ray 
beams. 

The axis of rotation is a precision commercial 
rotary table [11] having a radial concentricity and 
axial movement less than one micrometer. The 
table is driven by a stepping motor coupled to a 
threadless screw and a 48 cm long tangent arm with 
fine adjustment of the rotation obtained by 
piezoelectric transducers. The long first crystal is 



mounted in the center of the rotating table so that 
its diffracting planes are parallel to the axis of rota- 
tion to within a few seconds (of plane angle). 

Angles are measured by a polarization encoded 
Michelson interferometer having an angular sensi- 
tivity of a few X 10"^ seconds. The interferometer 
is shown schematically in Fig. 2 where the inset 
shows the over and under arrangement of the 
beams. The interferometer is illuminated by a com- 
mercial 633 nm HeNe laser [12] which emits two 
orthogonal linearly polarized frequencies sepa- 
rated by 1.8 MHz, The polarization sensitive beam 
splitter, b, divides the incoming light beam so that 
one polarization and frequency traverses one arm 
of the interferometer and the other polarization 
and frequency traverses the other arm. The beam 
splitter, steering elements (e, g, h), and compensat- 
ing element (f) are chosen and positioned so that 
the two interferometer arms have equal air and 
equal glass paths at zero angle. Elements c and d 
are 90* polarization rotators which rotate the plane 
of polarization of the outgoing and return beams 
respectively in order that both comer cubes are 
transverscd with the .same polarization orientation 
and that the output beams are directed away from 
the laser. The roof prisms, i and j, return the light 
to the corner cubes, k and I, doubling the angular 
sensitivity of the interferometer and keeping the 
output beams fixed in .space. The comer cubes are 
mounted on an arm which is rigidly attached to the 
rotating crystal table and retroreflect the light in 
each arm of the interferometer. 
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Fig. 2. Schematic diagram of Ihc angle inierfcrometer. a, c, g. 
h — beam sleering elemenls: b-polarizalion sensilive beam 
spliller; c. d— 90° polarization rotators; f— glass path compcri- 
saling plate: i. j — roof prisms; k, I— corner cubes. The insert 
ihows the passage of ihe beam through ihc comer cube— roof 
prism part of ihc ititcrfcrometer. 
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The interferometer output beam which includes 
the two orthogonally polarized frequencies is ana- 
lyzed by a 45' polarizer included in a commercial 
detector which provides a measurement difference 
frequency signal. The laser and its electronics also 
provides a reference difference frequem^ signal 
for the light emitted by the laser. For a stationary 
interferometer the two difference frequency signals 
are identical except for a phase difference 0, which 
is proportional to the difference in the optical path 
lengths of the two interferometer arms. Rotation of 
the interferometer arm causes a phase shift of the 
measurement difference frequency signal relative 
to the reference difference frequency signal. Con- 
versely, phase comparison of the measurement dif- 
ference frequency signal with the reference 
difference frequency signal permits measurement 
of the angular rotation. 

The 1.8 MHz difference frequency signals must 
be converted into signals which permit up-down 
counting for clockwise and counter clockwise 
rotation of the axis, measurement of fringe frac- 
tions, and servo control of the axis position. 
Although electronics to perform these functions is 
commercially available, no existing system could be 
conveniently coupled to our spectrometer. An elec- 
tronics module to perform the above functions has 
been developed and is briefly described in 
Appendix B, 

The diffraction profiles are recorded as x-ray in- 
tensity vs interferometer fringes and, thus, the an- 
gular offsets, ^are also measured in interferometer 
fringes. However, the angular offsets, p, needed in 
Eq. (2) must be in absolute angles. In order to con- 
vert interferometer fringes into absolute angles, the 
angle interferometer is calibrated using an optical 
polygon and a photoelectric autocollimator [4], A 
24-sided optical polygon (external angles —15°) is 
mounted on top of an indexing table on the rota- 
tion axis in place of the first crystal. The photoelec- 
tric autocollimator senses the directions normal to 
the faces of the optical polygon. The 24 external 
angles are measured in terms of interferometer 
fringes and the sum of the angles is constrained to 
be 360°. The equation relating angles and fringes is 
F~Ksm(6) where F is the fringe number, 8 is the 
corresponding angle of rotation, and K is the cali- 
bration constant. The measured value of K is 
5138551.7 at 22.5 "C. Near ^ = 0°, the value for the 
angular rotation/fringe is 0.0401 arcsec/fringe. The 
calibration constant is directly related to the sepa- 
ration of the corner cubes. Because the arm which 
defines the separation of the corner cubes is made 
oi]t of stainless steel, the calibration constant has a 



large temperature coefficient { — 16.7 ppm/°C). As 
long as the angular offsets are small and the tem- 
perature is well controlled, this large temperature 
coefficient can be tolerated. However, in the near 
future the spectrometer will be equipped with an 
invar corner cube arm. 

The zero angle of the interferometer is deter- 
mined by measuring one of the external angles of 
the polygon in a symmetric and asymmetric fash- 
ion. In the symmetric measurement, the polygon is 
positioned on the axis so that the faces forming the 
15" angle are measured at +7.5° and -7.5°. This 
symmetric measurement of the external angle is 
quite insensitive to the zero angle of the interfer- 
ometer. In the asymmetric measurement, the 
polygon is positioned so that the faces of the same 
angle are measured at ±14.5° and *0.5°. This 
as)Tnmetric measurement of the external angle is 
very sensitive to the zero angle of the interferome- 
ter. The zero angle is determined by requiring 
equality of the two measurements of the same an- 
gle. The uncertainty of the zero angle deter- 
mination is approximately 2 arcsec which makes a 
completely negligible contribution to the measure- 
ment of the small angular offsets near zero angle. 

The second crystals reside on a precision transla- 
tor equipped with a stepping motor and a linear 
encoder [13]. The slide has 15 cm of travel and a 
positioning accuracy of 0.01 mm. Although the 
slide has pitch and yaw of a few seconds in a few 
cm of travel, it is reproducible at the 0.01 second 
level and shows no short term drift. Each second 
crystal is mounted on a flexure hinge so its planes 
can be oriented to be roughly parallel to the planes 
of the first crystal. Fine adjustment is achieved by 
applying voltage to a piezoelectric "tipper." 

The diffracted x-rays are detected by two identi- 
cal Nal(Tl) detectors which pivot about the second 
crystal position. In front of each detector are three 
slits each 1,6 mmx 1.6 mm. The slits are arranged 
in the direction normal to the plane of dispersion 
with center of the middle opening on the plane of 
dispersion and the top and bottom openings 3.2 
mm above and below the plane of dispersion. Each 
.slit has a computer controlled shutter which per- 
mits radiation to pass through any one of the three 
openings. 

4. Crystal Preparation 

In order to make measurements with an uncer- 
tainty of 0.01 ppm, considerable care must be 
exercised in cutting, polishing, mounting, and 
aligning the crystal samples. The full width at half 
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maximum (FWHM) of nondispersive x-ray profiles 
obtained with Ag Ka radiation and the Si 440 re- 
flection is about 0.6 seconds. Since the angular sep- 
aration which corresponds to a lattice spacing 
difference of 1 x 10"" is —0.001 seconds, it appears 
that peak positions must be determined to 1/500 of 
the FWHM, This is a formidable task requiring 
very reproducible profiles. However, by carefully 
matching the thickness of the first and second crys- 
tals to within a few micrometers, nondispersive 
Laue-case rocking curves exhibit oscillatory fine 
structure which has a modulation period which is 
typically less than 1/10 of the FWHM of the profile. 
This fine structure is predicted by the dynamical 
theory of x-ray diffraction and has been used in 
high precision determinations of the structure fac- 
tor of Si [14, 15] and Ge [16] and in x-ray refractive 
index measurents [17]. We use this fine structure 
on the x-ray profiles as a sharp convenient refer- 
ence to measure the angular separation of the pro- 
files. 

The integrated reflectivity for the double crystal 
profiles oscillates as a function of crystal thickness. 
In Fig. 3, the integrated reflectivity vs crystal thick- 
ness is shown for the Si 440 diffraction of Ag Ka 
radiation. We have chosen a crystal thickness of 
0.455 mm as a compromise of intensity and practi- 
cal problems with thin crystals. Theoretical and 
experimental profiles for this thickness are shown 
in Fig. 4. 

The theoretical profile [18, 19] is obtained by 
adding the convolution of two intrinsic reflections 
for the two polarization states (<t and ir) 
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Fig. 4. Theoretical (lower) and experimental (upper) profiles 
for the 440 reflection of Ag Ka radiation using 0.455 mm thick 
Si crystals. 

where zl = the angle between the diffracting 
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Fig. 3. Integrated reflectivity of double crystal nondispersive 
profiles as a function of crystal thickness for the Si 440 reflection 
and Ag Kai radiation. 



ftj is the Bragg angle 

p. is the linear absorption coefficient 

t is the crystal thickness 

h is the extinction length 

y is the angular parameter 

C is the polarization factor 

A is the wavelength 

Xh is the complex electronic susceptibility of the 

crystal: }<t,^+ipCi,\ 

For low absorption where xn, > xu, it can be 

.shown that 

v^ ^ 1 - i (e/^o)/(-;3-cos0B) 
h =i7-Kcos&/(A/-J^|C|) 
r-ie) = [exp( - ^/cosft,)/(y^ + 1)] 
{[sin^[^ (V^+ 1)""] + sinA^[^e(y^+ l)-"^/2cosfl„]} 
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where e = (|x;Ml/ij£[!i!)|C[ 
Xhi, Xd* = imaginary parts of the Fourier coeffi- 
cients ?<t„ Xa 
K= volume of the unit eel! 
/'e = classical electron radius 
F= temperature modified structure factor 

When generating a theoretical profile, numerical 
values must be assigned to the following quantities 
ft., I, F, A, fc, y, and e. Conversely, when fitting an 
experimental profile with the theoretical descrip- 
tion, a quantity such as / or F can be varied until 
the difference between the experimental and theo- 
retical profiles is a minimum. 

In order to obtain strain free samples, the crys- 
tals are made with a 1.5 cm X 1,5 cm base on top of 
which is the thin 1.5 cm high wafer which is used 
for diffraction. Drawings of the long first crystal 
and the second ciystal are shown in Figs. 5 and 6. 
Typical lengths of the long and short crystals are 
6.3 cm and 1.2 cm, respectively. The long ct^^tal is 
balanced on a 1.5 nmi wde silicon rod which is 
centered on an optically poli,shed silicon cylinder. 
This assembly is waxed together with a low temper- 
ature optician's wax. The polished silicon c>'linder 
is attached to a ground cast iron base with epoxy 
resin. The cast iron base and the rotating table are 
joined by screw fasteners. After assembly the angu- 
lar offset between the normal to the polished sili- 
con cylinder and the crystal planes is measured 
using an x-ray spectrometer and an autocolllmator. 
The long crj'stal is mounted on a flexure hinge on 
the rotating table and is tipped an amount equal to 
the measured offset. In this way the crystal planes 
are made parallel to the axis of rotation with an 
estimated uncertainty of 2 seconds. 

The short crystal assembly consists of a ground 
cast iron base, a piezoelectric tipper, a silicon anvil 
(with a small hole for a thermistor), and the crystal. 
The base, PZT elements, and silicon anvil arc at- 
tached by conductive epoxy and the crystal is bal- 
anced on the raised portion of the silicon anvil (1^ 
mm wide) and attached by low temperature opti- 
cian's wax. This assembly is mounted on a flexure 
hinge on the translation table. Using x rays dif- 
fracted by the first crystal above and below the 
plane of dispersion, the planes of the second crys- 
tal are made parallel to those of the first crystal by 
the coarse mechanical adjustment and the fine 
PZT tipper. 

The crystals are cut with a diamond saw and 
then polished to the desired thickness using a 
chemically assisted mechanical polishing solution 
[20]. Typical polishing rates are 10 fim per h and 



50 iim to 70 fi,m are removed from each side t 
insure that no saw damage remains. Crystal thicl 
ness was measured using a coordinate measurin 
machine [21]. The reproducibility of the thicknes 
measuring procedure was 1 p.m-2 p.m. Polish in 
the long first crjstal so that the two areas used fo 
diffraction had equal thickness was the most diff: 
cult task. Because the edges of the crystals tend t 
be too thin, the areas used for diffraction were a 
least 2 mm from any crystal edge. After polishing 
all surfaces of the crystal except for the polishei 
surfaces were etched in HF/HNO3 to relieve strair 




FiB. S. Detailed drawings of the first ctyMal. a — cast iron, b- 
Si ilisk, and c— Si rod, and d-Si diffraction crystal. 




Fig. fi. IJciailed drawing of the second crystal, a — cast iron 
b — PZT, c — Si suppon, d — Si diffraction crysial. 
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5. Temperature Measurement and Con- 
trol 

Because the lattice parameter measurements are 
relative measurements, we required only accurate 
relative temperature values and not accurate abso- 
lute values. Since the expansion coefficient for sili- 
con is 2.56x10"*/K, relative temperature 
measurements accurate to 4 mK suffice for 1 x 10'* 
lattice parameter measurements. Six ultra-stable 
calibrated thermistor probes are used to measure 
the temperature of the crystals: two thermistors on 
the base of the long first crystal (one on each side) 
and one thermistor on the silicon anvil to which 
each of the four silicon crystals is attached. Vac- 
uum grease is used to insure good thermal contact 
between the silicon and the thermistors. The six 
thermistors and an ultra-stable calibrated standard 
resistor (2.5 kfl) [22] are connected in series and 
powered with a constant current source ( = lxlO~* 
A). A precision digital voltmeter is used to read the 
voltage drop across each thermistor and the stan- 
dard resistor. The temperatures are read every 4 s, 
so the temperature for a measurement time of 20 s 
is an average of five measurements. 

Before and after a critical measurement, the sLx 
thermistors are placed in a constant temperature 
bath in order to measure offsets be^veen the ther- 
mistors. The offsets are typically a few mK and are 
stable to less than 1 mK over times long compared 
to the few weeks needed for a lattice parameter 
measurement. 

The temperature of the laboratory is constant to 
within —0.1 K over several days. For the measure- 
ment reported here the spectrometer was isolated 
from direct air currents within the laboratory by ii 
cloth curtain. The curtain has subsequently been 
replaced by an insulated chamber. The crystals are 
covered with an aluminum thermal shield and two 
20 mm thick styrofoam thermal shields. The x-ray 
tubes, motors, and detectors are outside the 
shields. The temperature difference between the 
first and second crystals is typically 0.1 K and can 
he varied by changing the temperature of the water 
used to cool the x-ray tubes. The temperature dif- 
ferences between the second crystals is only a few 
mK. Variation of the temperature differences (the 
critical quantity for relative measurements) over a 
24 h period is typically less than 20 mK. 



6. Data Analysis 

The angular offsets, ^s and /3u, which are needed 
to calculate Ldid from Eq. (2) are the angular 



separation of profiles recorded along the path ( - ) 
and path ( + ). The profiles are measured by step- 
ping the axis through N discrete angles (typically 
N -WO) which are recorded as angle interferome- 
ter fringe numbers. For each angle the number of 
x-ray photons in each detector is counted for a 
fixed time (typically 10 s to 20 s). The time per scan 
is 40 min to 80 min. The fringe numbers for each 
profile are corrected in a point by point manner for 
the temperature difference between the first and 
second crystals. From the Bragg equation the de- 
pendence of the diffraction angle, &, on tempera- 
lure, T, is easily seen to be 

de/dr--«,tantf , 

where 00= linear coefficient of thermal expansion 
of the crystal. For silicon a»i = 2.56x lO'^/K. Let 8- 
and Oy be the angular settings of the first crystal 
for the path ( — ) and path ( + ), respectively, the 
rotation angle reference be the y axis, and clock- 
wise and counter clockwise rotations be positive 
and negative respectively. When the Bragg condi- 
tion is simultaneously satisfied at crystals 1 and 2 
and both crystals are at a reference temperature, 
Til, (sec Fig, la and Fig. 8), 



and 



0-(7-o) = 0b:(7'o)-M7'o) 



fl.(7'„)=-0B:(7'o)+0Hi(7'o). 



When crystal 1 is at temperature, Ti, and crystal 2 
is at temperature, T2, then these two equations be- 
come 



fl- (r„) = [M7"2) - ftjiCr,)] - ^ (7, - r,) = 



6e 



0-(TuT2)-^(T,-TO 



and 



e. (To)=-eB2iT2)+6si(Tt) + ^{T2-T,) = 



0ATuT2) + ^(T2-T,), 

where we have assumed the crystals are of the 
same material. 
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The angles 0- for the path ( - ) are corrected by 
-(\d/dT{T2-T,) while the angles fl+ for the path 
( + ) are corrected by -i-dd/dT (T^-Tt). Since the 
profiles are recorded in fringe numhcTs, these an- 
gular corrections need to be converted to fringes 
before adding. The angle measuring interferometer 
is positioned withiit a few minutes of the zero de- 
gree angle for all the measured offsets. Using the 
interferometer equation discussed above, the angle 
to fringe conversion factor of 0.0401 sec/fringe at 
zero degree angle is determined and used to con- 
vert the angle corrections due to lemperature dif- 
ferences into fringes. 

The angular offsets betv^'een the profiles on the 
temperature corrected fringe scale have been de- 
termined by two different methods. In the first 
method the profiles were fit with a iojrentzian 
function using a nonlinear least squares procedure 
in which the position, intensity, width, and back- 
ground are adjusted. The data with its pronounced 
wiggles and central spike i.s not well represented by 
the smooth Lorentzian function. However, by tak- 
ing the differences, /?=, between the recorded pro- 
files, /■,, and the smooth fitting functions, f.-., the 
wiggles of the two profiles are empha.<;tzcd. 

/?.(ft) = (/.(ft)-/-4ft)) 

The correlation function, C(A8), for the two sets 
of residuals is computed. 



i-i 

where the sum is over the N data points of the 
profile. 

The value of Ad for which C(A0) is a maximum 
is the angular separation (offset angle p) of the 
two profiles. The uncertainty in the measurement 
of the angular separation of the profiles is typically 
«6xlO-Vad, 

In the second method the dynamical diffraction 
function was fit to the recorded profiles. The crys- 
tal thicknesses were fixed at 0.455 mm and the 
structure factors were taken from Refs. [14 and 
15]. The only adjustable parameters in the fitting 
procedure were the position, the intensity and the 
background. The angular separation between tiie 
profiles is obtained as the difference between the 
fitted position of the two profiles. The measured 
angular separations obtained with the two methods 
agree within the measurement uncertainty. Be- 
cause the first method is computationally simpler 



and more easily adapted to a small computer, it 
was used to obtain all of the results presented 
below. 

A single lattice comparison run usually consists 
of 16 to 20 data scans, preceded and followed by 
alignment scans. The profiles are scanned by 
repeating (he following sequence: unknown crys- 
tal— cw rotation, unknown crystal— ccw rotation, 
standard crystal — cw rotation, and standard crys- 
tal— ccw rotation. After determining the angular 
separation for each .scan, the angular separations vs 
time of day for each crystal are fit with a variable 
order pol>'nomiaI (usually 3). The constant (in 
time) angular offset between these two curves is 
the angular offset between the unknown and stan- 
dard crystals for one data set. The scatter in the 
offsets is typically <6x 10""' rad provided the crys- 
tals are well aligned and the same areas of the crys- 
tals are used to diflraci the x-ray beam. 

7. Comparison of Two Silicon Samples 

Two Si samples were prepared as described in 
Sec. 4 from material acquired by PTB and NIST 
for absolute lattice parameter measurements. The 
NIST sample was supplied by Dow Chemical' and 
was a slab adjacent to the NIST x-ray/optical inter- 
ferometer. 'I"he PTB material was supplied by 
Wacker-Chemitronic and labeled WASO 17 by 
PTB. The relative difference between WASO 17 
and WASO 4.2 (the PTB x-ray/optical interferome- 
ter crystal) has been measured to be 2.5±1 x 10"** 
[23](WASO 17 -WASO 4.2). The first long crystal 
was prepared from material supplied by Monsanto. 

These two samples were chosen because the lat- 
tice spacing of the NIST sample has been mea- 
sured aksolutely and the lattice spacing of WASO 
17 was measured relative to the WASO 4,2 sample 
which was measured absolutely at the PTB, In the 
case of the NIST sample, further refinement of the 
absolute lattice parameter measurement is continu- 
ing [24, Z^l- Because the absolute lattice measure- 
ments have not yet achieved 0.01 ppm uncertainty, 
they do not provide as definitive a test of the Ad 
measurements as i,s desirable. The published re- 
sults of the absolute lattice spacing measurements 
are 



- Crorlfiin i:i.)mmercii»l i:qiiipmi;rt, instruiiitnts, iir materials arc 
iikntified in this papur lo specify adi:qual«ly the cxpcrinii;nt<)l 
procTdurt;. Such idemificalion doc^ not impl)' rCi:nminem1ation 
<ir endorsement by ihit National Instityte of .Standards and 
TccluKili:^', ntir dtjcs il imply that the materials or equipment 
id(;ntificd arc necessarily the hu-si available frir the purpose. 
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£/22d(WASO 4.2) at 22.5 "0 = 0.192015560 (12) nm 
(±6.2x10-*) 

d22o(NIST) at 22.5 "'0 = 0.19201554 (2) nm 
(±10.4x10-*) 

which lead to 
Arf/rf(NIST-WASO 4.2) = - 1 x 10-''± 1.2 x 10''. 

In order to develop confidence in the Ad mea- 
surements, we have measured the long term repro- 
ducibility by recording an extensive set of data over 
more than one half year during which the crystals 
were realigned on the spectrometer several times. 
In addition, measurements were taken using two 
different wavelengths (Ag Ka and Mo Ka x-ray 
radiation) and with the long first crystal in the two 
possible orientations (the 0° and the 180° orienta- 
tion). 

In Fig. 7 the measured values of Arf/rf are plot- 
ted vs time with different symbols for the different 
sources and different first crystal orientations. 
Note the long time span over which the Ag Ka, 0° 
orientation measurements were taken. Fn Table 1 
the average numerical values for the different 
sources and first crystal orientations are presented. 
From both the table and the figure, the magnitude 
of systematic effects related to the source and first 
crystal orientation is estimated to be —7x10"'. 
Table 2 provides a summaiy of uncenainty contri- 
butions resulting from systematic effects associated 
with 1) the wavelength and first crystal orientation, 
2) crystal temperature measurements, 3) crystal 
misalignments, and 4) periodic nonlinearity in the 
Michelson angle interferometer. 



Table I. Comparison of NIST and WASO 17 samples 



2 tJB 



•0-5 _ 









Tl«f (<"*») 



Rg. 7. Comparison of NIST and WASO 17 crystals: + AgKor, 
firsi ci>'stal 0°; D AgKa, first ciystal ISO"; O MoKn, first cn-stal 
0°; A MoKd, first cryslal 180°. 



Source 



First cryslal 
orientation 



M/d X 10" 
NIST-WASO 17 



AgKa 


0° 


0.508 + 0.379 


AgK« 


180° 


1.236 + 0.393 


MoKa 


0° 


1.050*0.262 


MoKa 


180° 


1.251 SOJIO 



T»bl« 2. Estimaled uncertainties from systematic effects 



Systematic effect 



Contribution to Arf tf x 10" 



A and 1st crystal orientation 
Temperature 
Crystal misalignment 
Interferometiy nonlinearity 



0.7 
0.3 
0.2 
OS 



The uncertainly in the measurement of the crys- 
tal temperature differences is less than 1 mK. This 
value was determined by periodically placing all of 
the thermistors in a constant temperature bath as 
explained in Sec. 5. A 1 mK uncertainty in the crys- 
lal temperature differences leads to an uncertainty 
of 3 X 10"' in Ad/d. The uncertainty introduced by 
crystal misalignment is discussed in detail in Ap- 
pendix A and is estimated to be less than 2x 10"'. 

Heterodyne Michelson interferometers are prone 
to sub-periodic nonlinearity resulting from imper- 
fect separation of the two frequencies by the pri- 
mary beam splitter [26, 27]. The periodic 
nonlinearity can be evaluated by pressure scanning 
the interferometer and recording the fringe ad- 
vance vs pressure increase. The maximum ampli- 
tude of the nonlinearity is estimated to be 4 x 10"" 
rad. This means that the angular separation of two 
points on a diffraction profile which are separated 
by 0.5 interferometer fringe (4x 10"' rad) might be 
in error by 8 X 10-' rad. However, because the pro- 
files are typically 25 interferometer fringes wide, 
the influence of the interferometer periodic error 
on the profile peak position is significantly reduced 
by averaging. In addition, the phase relation of the 
x-ray profile to the angle interferometer fringes 
changes so that uncertainty due to periodic nonlin- 
earity becomes part of the statistical uncertainty. 
We estimate that the uncertainty contribution to 
Ad/d due to the systematic effect associated with 
the periodic nonlinearity of the angle interferome- 
ters is not larger than 5 x 10"'. The angle interfer- 
ometer is being modified in order to reduce the 
periodic nonlinearity by approximately a factor of 5. 
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A final value for ttie NIST-WASO 17 compari 
son was obtained by considering each of the entries 
in Table 1 to be an independent measurement of 
equal weight. The statistical uncertainty was com- 
bined with the uncertainties from systematic effects 
in Table 2 to obtain a value for Ad/d (NIST- 
WASO 17) = (1.037± I.O)x lO-l By combining this 
value with the WASO 17-WASO 4.2 difference 
noted above, one obtains 

Adfd (NIST-WASO 4.2) = {3 .5 ± 1.4) X 10". 

This value is only slightly outside the 1 tr uncer- 
tainty of the absolute measurements. In addition, it 
should be remembered that the NIST absolute 
value is a preliminary result. 



8. Conclusions 

A lattice comparison facility has been estab- 
lished at NIST which is capable of measuring ciys- 
tal lattice spacing difference with an uncertainty 
<0.01 ppm. By using crystals of equal thickness, 
the recorded profiles exhibit fine structure which 
permits more precise measurement of the small an- 
gular offsets between profiles. The spectrometer is 
designed to permit easy interchange of crystal sam- 
ples and the comparison of four samples in one 
setup. Procedures are provided for precise align- 
ment of the crystals so that alignment errors con- 
tribute < 2x 10"'. A comparison of samples whose 
lattice spacings have been absolutely measured is 



consistent with the absolute measurements. In tf 
near future other crystal .samples (including Gt 
destined for x- and gamma-ray diffraction will t 
compared using the spectrometer and the tecl 
nique described here. 

9. Appendix A. Derivation of th 
Equation for Calculation of Ad 

If the lattice planes in cr>'Stal 1 and 2 are a 
perfectly parallel with the axis of rotation and th 
x-ray beams lie in a plane normal to the axis c 
rotation, then the difference in lattice spacing bt 
tween crystal 1 and 2 can be derived from Bragg 
law (see Sec, 2). In reality, however, the alignmer 
of the crystals and/or the x-ray beams is never pei 
feet. In this appendix, equations for the calculatio 
of Ad which include misalignment terms ar 
derN'ed. The x-ray beams are described b 
wavevectors, Jt, and the crystals are described b 
reciprocal vectors, t. The derivation is based o: 
the fact that the angular positions, 0+ and 6^, c 
the symmetric rocking curve profiles, measurei 
along path ( + ) and path (-) in the two crysta 
spectrometer, are equal to the two particular rota 
tion angles 0, for which the Laue-conditions fo 
path { + ) and path ( — ), are simultaneously satis 
fled at crystal 1 and 2 (.see Figs. 8 and 9). 

Lar.son [28] has made a similar geometrical anal 
ysis of a lattice parameter comparison using Bragj 
geometry. We first derive the equation for Ad/d fo 
the ideal case involving perfect crystals and no mis 
alignment of the crystals or the x-ray beams. Tht 



'..+ 



T?^ 



-*^^^. 




Fir. S. Geometry of the lattice parameter measurement with no crystal mis- 
align tncnt. 
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Fig. 9. Geometry of the lattice parameter measurement including crystal 
misalignment. 



more experienced reader may choose to skip over 
this section. Next, we derive the equations for the 
nontdeal case and apply them to "alignment scans" 
and "data scans." 

9.1 Ideal Case 

First consider the ideal case of no misalignment 
of the x-ray beams or tilting of the crystals. The 
geometry of the measurement is then entirely con- 
fined to the xy plane and is as shown in Fig. 8. The 
coordinate system is such that the z-axis is parallel 
to the e-axis (i.e., the axis of the rotating table), 
thejc-axis is parallel to the lattice planes of interest 
in crystal 2, the y-axis is normal to the lattice 
planes in crystal 2, and the normal to the lattice 
planes in crystal 1 lies in thexy plane. The rotation 
angle reference is taken to be ihey-axis and the 
fl's, the angles of rotation of the first crystal, are 
positive for clockwise rotation and negative for 
counterclock\\ise rotation. 

The x-ray beams are described by their wave vec- 
tors: Jt., Jt^ for the incident rays; kL, k^ for the 
singly diffracted rays; and A", kl for the double 
diffracted rays. Using the angles shown in Fig. 8, 
the doubly diffracted rays can be written as 

it-=27rM(cosa-, sina', 0), 

Jt- =2Tr/A(cosa+, sina", 0). 

The crystals are described by their reciprocal lat- 
tice vectors Ti (0) for crystal 1 and n for crystal 2 
which can be written as 



27 



Ti(fi) = ^ (sine, cosfl, 0). 
7^=^ (0.1,0). 



The planes at the two ends of crystal 1 are assumed 
to be parallel and to have the same lattice spacing. 
When the Laue condition is satisfied at crystal 1 

k'.'-k-=-ni6.)andki-k.=T,(6^). (3) 

Similarly when the L^ue condition is satisfied at 

crystal 2 



kl -Jt: = n and kl -kl = -n. 



(4) 



Since the wave vectors have equal magnitude, 
2t7/A, Eqs. (3 and 4) imply 



Ti(0O ■ (7^"A:")^|Ti(e-)P and r,(0O - (7^+*:) = 
||T.(MP. (5) 

By substituting the reciprocal lattice vectors and 
the wave vector with a-=a'^. = fe = Bragg angle at 
the crystal 2 into Eq, (5) and using the Bragg equa- 
tion, A =2^2 sin 02, to eliminate the wavelength, two 
equations for the lattice spacing difference are 
obtained. 
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^^-7^= -sine- cot 0Ha-(l -cos0-)^ 
-^ — *-= -sinS+cotfe;-(l -cos0-) = 



e^COtflaj-r^^ 



where terms up to second order in 6^ have been 
retained. Tliese two equations can be added to ob- 
tain 

^ = _if^Z^ coteB,-|(0--0.)=, (6) 

where we have used the fact that flf = — 8-, 

For comparisons in which the first and second 
crystal are made out of tlie same material 
[(dz—di)ldi<iW~^\ the last term is negligible as 
the following inequality shows 



10 



-10- 



d2~d^ 



d. 



'>^ 



d^-dA^ 



tan"0B2 = 



i (e. -6^+)^if eB2<54°. 

If crystal 2 is alternately the standard and unknown 
crystal, the Eq. (6) becomes 



d,~di §, 

— S5S *-- 

di 2 



cotftjs where ^= fl-s— 5+, 



^■1 - tf I . Ai 



cotflau where pj= ff- „ - 8+a 



dt 2 

These two equations immediately lead to Eq. (2). 



with respect to the plane normal to the rotatioi 
axis introduce correction terms to the measurei 
lattice spacing difference. In the analysis whicl 
takes into account misalignments and crystal im 
perfections, the geometry of Fig. 9 is used. 

The coordinate system is identical to that use( 
for the ideal case, i.e., the z axis is parallel to th( 
axis of rotation and the at axis parallel to the lattic( 
planes in crystal 2. However, the wave and recipro 
cal lattice vectors no longer lie in the^y plane. Thi 
wave vector of the double diffracted x-rays can b< 
put in the following form: 



Jt± =— (cosa± COSSr, 



^sinar, costti sinSi), (8^ 



where S- is the angle between thejiy plane and the 
projection of the wavevector Ar= on the xz plane. 



The reciprocal lattice vectors a,ssociated with path 
( + ) and path ( — ) in crystal 1. t+ and t-, are nc 
longer parallel and of equal magnitude. The twc 
vectors 7^ can be described by the six parameten 
d±, p~ and 0~ as follows: 

2Tr 
Ti =— (sin(6-<^-) coSjr>±, 



cos(0 — ^i) cospi, sinp=), 



(9) 



where rf* is the lattice spacing in the diffracting re- 
gion of cry.stal I corresponding to path ( i ), p^ is 
the angle between t. and the Jty-plane, and <^^ is 
the particular angle of rotation 6 for which ti is 
parallel to the jz plane. For a perfect first crystal, 
d+=d~, p+=p-, and 0+=4'_. In crystal 2 the re- 
ciprocal lattice vector, n, is the same for both paths 
and is given by 



277 



T!=-j- (0, cosp, sinp). 



(10) 



d, ^ 1 



colons 



(7) 



9 J. Non-Ideal Case 

Misalignment of the planes of crystals I and 2 
with respect to the rotation axis, imperfections in 
the first crystal (i.e., lattice spacing gradient and 
nonparallelism of the planes at the two ends of the 
crystal), and the misalignment of the x-ray beams 



where di is the lattice spacing in the diffracting re- 
gion of crystal 2, and p is the angle between -n and 
.vy-plane. Note that n is in theyz plane. 

For path (±), the relations between the wave 
vectors k~ of the incident rays, kL of the single dif- 
fracted rays, and kZ of the double diffracted rays 
are given by the Laue conditions: 



kL~k± = ±Ti (crj'stal 1) 



k'^—kl—-i-Ti (crystal 2) 



(11) 



(12) 
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Since the wave vectors have equal magnitude, 
lirl\, Eqs. (11) and (12) imply that 



UK _ 1 I 12 



r._-{-n±kV)=^r.\\ 



(13) 
(14) 



When the Laue condition is satisfied at crystal 2, 
the angles, a±, can be expressed as art = &B2+£± 
where €± are the small angular corrections to the 
second crystal Bragg angle because of crystal imper- 
fections and misalignments. Thus Eq. (8) becomes 



*t=~ (cos(&2 + e^) 



cosSi 



:sin(fe + ci), cos(0B2+ei) sin5±). 



(15) 



In order to establish a relation between the lat- 
tice spacing d of crystal 2 and the rotation angles 
= 0±, for which the Laue-conditions correspond- 
ing to path ( ± ) are simultaneously satisfied at crys- 
tal 1 and crystal 2, Eqs. (9), (10), and (15) are 
expanded to second order in the angles 6± — 4>±,p±, 
p and Sx and to first order in the correction angle 



r.=^je.-<f>.. 



i-lie.-<f>..y-lpl,p.)iu) 



T. = f (0,1 4 p=,p) 



(17) 



*:£=j((l-|5|)cot6>B2-ei, ^l + eiCOtfez, 

5.cote„,). (18) 



In the last equation, the wavelength A has been 
eliminated using A =2<i2sin0B2- By insertion of Eqs. 
(17) and (18) into Eq. (13), an expression for the 
correction angles e± is easily obtained: 



1 



€±=r p tan6'B±p5± 



(19) 



Now by inserting Eqs. (16), (17), (18), and (19) into 
Eq. (14) the following equation can be derived: 



d2-d^ 
dx 



= ±(^--01) COtfe + Ci , 



(20) 



where e* is the second order correction term 



e. = -^(0.-<^.)^-i(p.-p)(p.-3p) 



1 
2 



±5i (px-p) coXOez 



(21) 



The angle S± in Eq. (21) describes the inclination of 
the double diffracted beam relative to the xy -plane, 
(see Fig. 9). If crystal 1 and crystal 2 are perfectly 
parallel, the wavevectors Jt±, ki, and kx are copla- 
nar and the inclinations S± are determined simply 
by the positions of the source and detector slits. 
When the two crystals are not perfectly parallel, d± 
are also dependent on the misalignment angles p± 
and p. 

Let €, €', and C be projections along thex-axis of 
the x-ray paths from the source slit to crystal 1, from 
crystal 1 to crystal 2, and from crystal 2 to the detec- 
tor slit, respectively. The distances are shown in Fig. 
1 and for our spectrometer are: 

€ =(25.4 cos^ + 8.25) cm 

e' =8.25 cm 

r = (28.5 cos 0b) cm 

If AZi denotes the difference in Z- coordinate be- 
tween the detector slit and source slit of path (±), 
then the inclination 4+ relative to the xy-plane of 
the straight line joining the two slits is given by: 



Ax^ 



AZ. 



^ e+e'+r- 



(22) 



The change in Z -coordinate as one moves from the 
source slit to the detector slit along path ( ± ) can be 
expressed by combining the Z-components of the 
wave vectors kx, ki, and JfcC, normalized to unity 
and the distances traveled: 

AZ± =^ [(k^).e + ikl)J' + (Jt;).r] 



1 
cosfe ■ 



(23) 
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From Eqs. (11), (12), and (15) and using Eqs. (9) 
and (10), the Z-componenls of the wave vectors 
are found to second order in p*, p, and 5±: 



^:;-(it=).=5:. cosft)±2(p-pi)sin£?B (24) 



-^ {k'±). = 5= cosAb ± 2p sinfle (25) 



-^(kl)._ = S^cosfk. (26) 

Air 



By insertion of Eqs. (22), (24), (25), and (26) into 
Eq. (23) a relation between 5. and A 5, is obtained: 



Using this expression, the unknown angles 8± can 
be eliminated from Eq. (21): 

e±=-|(0±-<&O' + 2(p.-p)(a,Pi+6,p)± 

il±(pi'p) cot(fe (27) 



where 



°^=e+t'^r-r^-'^ 



This term is two orders of magnitude smaller than 
the intended uncertainty (a few parts in lO**) in the 
lattice spacing comparison and is therefore ne- 
glected in the following. The remaining two terms 
in Eq. (27) are both zero, if the lattice planes in 
crystal 1 and 2 are parallel, i.e., p+— p = and 
p.-p=0. The aim of the alignment scan is to 
measure p+—p and p- — p and then to adjust the 
tip of crystal 2 so as to make p*—p and p,~p as 
close to zero as possible. 

9.2.1 The Alignment Scan In the alignment 
scan, only one path, path (+ ) or path ( — ), is used. 
Two rocking curves are measured simultaneously: 
one curve using the upper (up) position and one 
using the lower (lo) position of the detector slit. 
For the path (-), the observed difference A0- in 
the centers of the two rocking curves 
(A0- = 0.„p-0-io), the change in the inclination 
{A-uf — A-\a) when the detector slit is moved from 
the upper to the lower position, and the difference 
in misalignment of the two crystals, p-—p, are re- 
lated. The following equation is derived directly 
from Eqs. (20) and (27): 



P--p = 



-A9- 



A-up— 4-lo 



(28) 



Similarly, for path ( + ) it follows that 



P*-P = 



Ae* 



A+^— 4n 



(29) 



fct = 



r 1 

-4 = 0.15 



^-i-r-f-r 4 



Even though the distances t and t depend on the 
Bragg angle fc, the numerical values for a\ and b] 
are essentially identical for the two x-ray sources 
(fib = 16.9" for Ag Kaj and ft, = 21 .6= for Mo Koi). 
When crystal 1 and 2 are both of the same mate- 
rial, an upper limit of the first term in Eq. (27) is 
easily estimated: 



1 



\d~d.. 



2i^^--^^f<\-jr\ < 



10- 



10" 



k-dr 



: 1 

=■2 



—^ — tan'flB=2(^±~'^*) ' 



if 0B<54'. 



(A-^—A-io) and (A+up-4+io) can be estimated 
from mechanical measurements of the spectrome- 
ter and are ~9 X 10"^ rad. Due to bending of crys- 
tal 1, p+ and p- are not equal. Consequently, the 
sample crystal can only be aligned for one of the 
two paths. If the tilt of crystal 2 is adjusted so that 
\A8-\ =2x10'^ rad, then tj-pical values for \Ae+\ 
are 1 x 10"^ rad. These values imply that 

|p--p|=2xl0-'' 

|p+-p] = lx10-^ . 

The magnitude of the individual p's can also be 
controlled and estimated. By rotating the 1st crys- 
tal by 180°, the tilt of the first crystal can be ad- 
justed so that p+ and p- are < 1 xlO"' rad. 

9.2.2 The Data-Taking Scan In the data- 
taking scan, two rocking curves are simultaneously 
measured with the detector slits fixed in the middle 
position: one curve using path ( + ), and one using 
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path ( - ). We define the average lattice spacing, rfi, 
in crystal 1 as 

"^'-dTTdZ- 

and average Eqs. (20) and (27) over path ( + ) and 
path ( - ) to obtain 

^^^ = i (fl^ - e_ + r}u) cotflB +e, (30) 

where 
« ={(.P*-p)io\P*+bip) + (p--p)iaip- +bip)}+ 

~{A.{p.-p)-A.{p-~p)}coteB (31) 



and 



■i}o=4>- —4>*- 



(32) 



The angle tjo describes a bending of crystal 1 around 
the e-axis and is independent of the particular sam- 
ple used as crystal 2. Obviously, the appearance of 
the unknown angle 170 in Eq. (30) exclude the pos- 
sibility of using crystal 1 as an absolute reference 
crystal. 

The two crystals being compared are placed on 
the translation table and are used alternately as 
crystal 2. One of the crystals serves as the standard 
(s) and one as the unknown (u). The difference in 
lattice spacing. Ad =di — du for the two sample crys- 
tals can be derived from Eqs. (30) and (31). 

^=_(&zA)cotOH + Ae (33) 

where: 

Ac = [a; (/) * + p- ) + /Ji (ps -I- pu) - 2 

(4+-4-)cotft,](p,-A.) = 

[a2{p.+p-) + 2h2p„-^(^^-A.)cotdB] 

(ft - Pu) + 62 (ft - Pg)^ (34) 

where flj=b,-a,= -0.08 and b2= -2ii= -0.30. 

In order to make Ae small, the tilts of the crystals 
are adjusted so that p's are < 1 x 10"^ (see Sec. 



9.2.1). Care is taken to center the source and detec- 
tor slits on the plane of dispersion in order to make 
A+ — A~ small, but mechanical measurements sug- 
gest that the magnitude of the 4*s is ^ 3 x 10"*. The 
correction term can be estimated by measuring 
(ft— Ai) as a function of (ah— Ai) where p^ is kept 
fixed and p, is varied. For pt-pu S5xl0"^ Eqs. 
(33) and (34) can be approximated by 

^=-(^)cotfe + C(p.-p.) 

where 



C^^szip* +p-) + 2b2pu-2 (A*-A-)coteB. 

A plot of (pi - ft) vs (p, - /)u) is shown in Fig. 10 and 
has an approximate slope of — 3.6x 10"* in reason- 
able agreement with the estimates of the A's from 
the mechanical measurements and the estimates of 
the individual p's. 



K 



-C.B 



\ 



/.', ■■'].';. '.s !.3 ■l.^ 7-5 

Tn'^^ul X 10^ ' radians 

Fig. 10. Measured data profile offsets, ft -ft,, as a function of 
the measured crystal tilts (ft — ft,). Sec text for more explanation. 

Before and after data scans are recorded using a 
standard and unknown crystal, alignment scans are 
recorded for the standard and unknown crystal also. 
The quantity p, — pu is determined from the align- 
ment scans and is maintained at a value <4xl0"* 
rad. Thus the correction term Ae <2 x 10"'. Correc- 
tion of the data for measured values of Ae does not 
reduce the statistical spread of repeated measure- 
ments. We have thus chosen to increase the relative 
uncertainty of the final results by 2x 10"', but not 
to correct individual measurements. Although the 
correction term is small, it could be reduced by an 
order of magnitude by making the source and 
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detector slits adjustable about the plane of disper- 
sion. Modifications to the spectrometer to permit 
this adjustment are in progress. 

10. Appendix B. Electronics Associated 
with the Interferometer 

The layout of the Michelson angle inierferome- 
lers is shown in Fig. 2 and is briefly described in 
Sec, 3. In this appendix the electronics used to 
mea.'nire the whole and fractional fringes and to 
servo-control the rotation axes are described. A 
block diagram of the electronics is provided in Fig. 
11. 

The 1.8 MHz TTL level signals from the laser 
and the detector are first converted to low 
impedance signals by buffer amplifiers. The laser 
beat signal is then passed through a 2-way 90° split- 
ter with 0° and 90° split yielding sin (ajl) and cos 
(wf), respectively. The detector beat signal is 
passed through a 2-way 0° splitter yielding two sin 
(wf + <^) signals. The laser'.*; cos (wt) and one of 
the detector's sin ((ot + <f>) arc combined in a mixer 
which gives signals at the sum and difference fre- 
quencies. The difference frequency signal is sin ^ 
and the sum frequency signal is sin (2w( +<^) which 
is attenuated by a low pass filter. Meanwhile the 
sin (lift) from the laser is combined with the other 
sin (to/ + <^) of the detector in another mbter which 
after filtering gives cos <t». The sin ^ and cos i> sig- 
nals are next amplified to restore their signal 
strengths which have deteriorated due to losses in 
the splitters, mixers, and filters. 

Ninety degree changes in ({) which correspond to 
~0.04 second rotation of the rotary table are de- 
tected by sending the sin ^ and cos ^ signals to an 
up/down counter. Smaller changes in <J) are mea- 
sured by comparing <{> to a standard angular signal 



(sin and cos d) obtained as the output of a digital 
vector generator. A vector generator with a few 
0.1° sensitivity provides a total angular sensitivity of 
a few X 10"* seconds with the vector generator ad- 
dress being directly related to the fringe fraction. 
By taking the difference of the product of sin ^ and 
cos B and the product of cos<^ and sin 6 the signal 
sin(<^ - 6) is obtained. This is an appropriate error 
signal for locking the angle of the rotary table via a 
piezoelectric transducer to the position where 
4> = d. 



11. References 

[1] p. Becker, K. Dorenwendt, G. Etwling, R. Laucr, W. Lucas, 
R. Probst, H.-J. Radenmachcr, G. Reim, P. Seyfried, and 
H. Siegcrt. Absolute Measurcmeni of the (220) Lattice 
Plaue Spacing in a Silicon Crystal, Phys. Rev, Lett. 46, 
1540-154.1 (1981). 

[2J R. D. Dcslaiies and A. Henins, X-Ray to Visible Wave- 
length RatkK, Phys. Rev. l.ctt. 31, 972-975 (1973). 

f3] R. D. Dcslattcs, A. Renins, H. A, Bowman, R. M. 
Schoonovcr, C. L. Carroll, \. L. Barnes, L. A. Machlan, L. J. 
Moore, and W. R. Shields, Determination of the Avogadro 
Constant, Phys. Rev. Lett. 33. 463 (1974); R. D, Dcslattes, 
A. Hcnin.l. R. M. Schoonovcr, C. L. Carroll, and H. A. Bow- 
man, Avogadro Con slant -Correct ions to an Earlier Report, 
Phys. Rev. Utt. Ciimm. 36. Byit-900 (1976). 

[41 R- D. Deslaitcs, E_ G. Kcsslcr, W. C. Saudcr, and 
A. Hcnins, Remeasurcnncni of 7-Ray Reference Lines, 
Ann. Ph>-s. 129, 378 (1980)- 

[5] G. Basile, A. Bcrgamin, G. Cavagncro, G, Mana, E. Vit- 
tone, and G. Znsi, Silicon Lattice Constant: Limits in 
IMGC X-Ray/Optical Inlerfcromctry, IEEE Trans. In- 
stram. Mca.";. 40, 98-102 (lyVl). 

[6] M. Tanaka, K, Nakayama, and K. Kuroda, Experiment on 
the Ahsolutc Measurement of ii Silicon Lattice Spacing at 
the NRLM. IEEE Trans, Instrum. Meas. 3», 206-209 
(1989). 

[7] M. Hart, High Precision Lattice Parameter Measurements 
by Multiple Bragg Reflexion Diffractomctry, Proc. Roy. 
Soc. Ltmdon Ser. A 309, 281-296 (1969). 



S»^ , »»d3^i * ^i 



ANta-E 
Ih'TER- 
rCPOMETER 


■ - 


HECEiVER 


TBl»HI 


OF 
SPU1 ItH 
















1,8 MMl 


SPLITTEH 




LASER 




sitKvll 



^ 



.1 ._ 



_r 




nl'V-ft) 



COS'^-C'OTtOri t . f^| 



PULSE 
GEHEJIiTOB 


UP 
~ DOWl 


UPt>OWN 

ctumltlR 



■J VFCIOB 
^ j C*-?*E BATOR 



Fig. 11. Biiitk diagram nf cletlriirics used to measure whole and fractional fringes. 

17 



Volume 99, Number 1 , January-Fcbruaiy 1994 

Journal of Research of the National rnstitute of Standards and Technology 



[8] M. Ando. D. Batlcy, and M. Hart, A Simple Bra^- 
Spacins Comparator, Acta Oryst. A 34. 4S4-4S9 (I97S). 

|9] P. Becker, P, Seyfried, and H. Siegert, The Lattice 
Parameter of Highly Pure Silicon Single Ciy!>1al£, Z. Phvs. 
B 48, 17-21 (1982). 

[10] D. Haiu«rmann and M. Hart, A Fast High* Accuracy 
Lattice-Parameter Comparator, J. Appl, Cryst, 23, 63-69 
(1990). 

[11] The rotary table ix an "Ultraron Rotary Table," Model 
No, UR-6 supplied by A. A. Gage Inc., Detroit, Michigan. 

[12] The la^er is a Model SSOIA loser transducer, supplied by 
Hewlett Packard, Santa Clara, California. 

[13] The slide is a Model 7-6 Aitoride Tabic supplied by 
Anorad Corporation, HauppaUge, New Vorfc. 

[14] M. Deuisch and M. Han, Electronic Charge DistributkMi 
In Silicon, Phys. Rev. B 31. 3S4&-3858 (1985). 

(IS) R. Teworte and U. Bonsc, High-Precision Determination 
of Structure Faaors F» of Silicon, Pbys. Rev. B 29, 2102- 
2108 (1984). 

[16] M. Deut^ch, M. Hart, and 5. Cummjngs, Hi^-Accuracy 
Structure-Factor Measurements in Germanium, Phy. Rev. 
B 42, 1248-1253 (1990). 

[17] M. Deutsch and M. Hart, X-Ray Refractive- Index Mea- 
surement in Silicon and Lithium Fluoride, Phy. Rev. B 30, 
640-642 (1984). 

[18] U. Bonsc, W. Gracff, R. Tcworic, and H. Rauch, Oscilla- 
tory Structure of Laue Case Rocking Curves, Phy, Stat. 
Sol. 43. 487-492 (1977). 

[19] U. Bonse, and W, Graclf, X-Ray and Neutron Interfcr- 
ometry, in X-Ray Optics. H. J, QucLs.<ier. ed., Springer, 
Berlin (1977), pp. 9J-143. 

[20] The grinding solution is CAB-O-SPERSE SC-3010 sup- 
plied by Cabot Corpora titHi, Tusct^. Illinois, 

[21] Wc arc indebted to B. Borchert. R. Vcale, and W. 
Gallagher of the NIST Automated Manufacturing and 
Research Divtston for making cciordinatc measuring ma- 
chines available for these measurements. 

[22] Wc are indebted to N. Belecki and R. Green of the NiST 
Electricity Division for calibrating the stiindard resistor, 

[23] D. Windisch and P. Becker, Silicon Lattice Parameters as 
an Absolute Scale of Length for High Precision Measure- 
ments of Fundamental Constants, Phys, Stat. Sol. (a) IIS, 
379-388 (1990). 

[24] R. D. Dcslattcs, X-Ray Interferomelry and y-Ray Wave- 
lengths, in The Art of Measurement, B. Kramer, ed., 
Wcinheim, Germany; VCH Verlogsgescllschaft (1988), 
pp. 193-2(}8. 

[25] R. D. Deslaties and E, G. Kcsslcr, Jr., Status of a Silicon 
Lattice Measurement and Disseminatioo Exercise. IEEE 
Trans. Instnim. Meas. 40, 92^97 (1991), 

[26] C M. Sutton, Nonlinearityin Length Measurement Using 
Hetcrodyrw Laser Michelson Inteiferometry, J. Phys. E: 
Sci. Instram, 20, 1290-1292 (1987). 

[27] N, Bobnrff, Residual Errors in Laser Inierferometry from 
Air Turbulence and Nonlinearity, Appl, Opt, 26. 2676- 
26S2 (1987). 

[23] B. C Larson, High-Precis bn Measurements of Lattice 
Parameter Chnngcs in Neutron-Irradiated Copper, J. 
Appl. Phys. 45, 514-51K (1974), 



About the authors: Ernest G. Kessler and Albert 
Henins are physicists in the Quantum Metrology Divi- 
sion of the NfST Physics Laboratory, Richard D. 
Deslaties is a Senior NfST Research Fellow and Chief 
of the Quantum Metrology Division of the NJST 
Physics Laboratory. Lars Nielsen is employed by the 
Danish Institute of Fundamental Metrology and par- 
ticipated in this work when he was a guest scientist in 
tfte Quantum Metrology Division, Mohammad Arif, a 
physicist, was in the Quantum Metrology Division 
when this work wad done and is now in the Ionizing 
Radiation Division of the NIST Ph)'sics Laboratory. 
The National Institute of Standards and Technolo^ 
is an a@zncy of the Technology Administration, U.S. 
Department of Commerce. 



i$ 



